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Abstract. The spherically symmetric perturbations in the spatially flat Friedman 
models are considered. It is assumed that the Fricdmannian density and pressure 
are related through a linear equation of state. The perturbation is joined smoothly 
with an unperturbed Friedmann's background at the sound horizon of perturbation. 
Such junction is in accordance with the "birth" of a local perturbation as a result 
of the redistribution of matter. The solution of the Einstein's equations is obtained 
in linear approximation on a Friedmann's background near the the sound horizon of 
perturbation. 



PACS numbers: 04.25.Nx, 04.40.-b, 04.40.Nr 



Local spherically symmetric perturbations of spatially flat Friedmann models 



2 



1. Introduction 

The description of metric perturbations in relativistic cosmology is significance in 
the theory of creation of the universe large-scale structure. Spherically symmetric 
perturbations are an important and interesting class of that perturbations. 

The study of small perturbations on the cosmological background has a long history. 
The first relativistic treatment was given by Lifshitz T and developed by himself 
and Khalatnikov |2j. The gauge-invariant formalism of perturbation was given in the 
fundamental paper of Bardeen [3] and was developed in the works 0j. 

Spherically symmetry perturbations on the cosmological background have been 
studied by a lot of authors El EI- As it was noted in [7] the spherically symmetric 
perturbation in fluid with nonzero pressure may be represent as outgoing and ingoing 
waves (travelling from and into the center of the configuration with the velocity of 
sound). The ingoing wave amplitude increases as soon as the wave approaches to 
the center of the configuration. Therefore the magnitude of spherically symmetric 
perturbation near the center of configuration must increase more quickly than the 
amplitude of Lifshitz's plane-symmetric wave perturbations pQ. Moreover the local 
perturbation must be smoothly jointed with an unperturbed Friedmann's background 
at the sound horizon of perturbation. 

Our main difference from the paper by Ignat'ev and Popov (Zj is that we consider a 
cosmological fluid with a linear equation of state e = 7P, where e is density, p is pressure 
of fluid but leave arbitrary a constant 7 (7 > 0). 

In Sec. II we write out the rigorous equations describing the perfect fluid in 
spherically symmetric space-times and give the characterization a background space- 
time. In Sec. Ill we fix the gauge and set out the basic equations governing the 
spherically symmetric perturbations. The particular solutions of these equations are 
given in Sec. IV. In Sec. V we impose the boundary conditions on the perturbation and 
obtain the solution for perturbation near the boundary. Finally, the main conclusions 
raised in this paper are summarized in Sec. VI. 

2. Perfect fluid in spherically symmetric space-time. Rigorous equations 
and background space-time 

With the gravitational field being spherically symmetric, coordinates may be chosen so 
that the metric takes the form 



where v, \i and A are functions of the time coordinate r\ and a radial space coordinate 
r. As we are interested in perfect-fluid distribution, the energy-momentum tensor is 



ds 2 = -e v dr] 2 + e x dr 2 + eV 



(d9 2 + sm 2 9 dip 2 ) 



(1) 




(2) 
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where e is density, p is pressure and ■u M (-u r? , u r , 0, 0) is four-velocity of fluid. Take into 
account the equation of state 



e = 7p 

one can obtain the following nontrivial Einstein equations (c — G—1) 



(3) 
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where v = u r e x ^ 2 is the frame projection of radial velocity, a dot and a prime denote 
partial derivatives with respect to rj and r, respectively 

In isotropic spherical coordinates the metric of the background space-time can be 
written as 



ds > = a 2^ dr] 2 + dr 2 + r 2 ( d0 2 + gin 2 Q ^2) } ? 

where the function 0(77) satisfy the equation 



(8) 



2 7 - = (7 - 3) - (9) 

a cr 

The energy density e , the pressure p and the frame projection of radial velocity of the 
background fluid v are 

■ 2 



£0 



3 a 
8na 4 



, £0 = IPo, v = 0. 



(10) 



3. Gauge transformation and perturbed equations 



We shall only consider spherically symmetric perturbations of spatially flat Friedmann 
models, i.e. 



v = In (a 2 ) + Su, X = In (a 2 ) + 5X, p = In (a 2 ) + <5/i, 
£ = e + Se, p = Po + 5p, 



(11) 
(12) 
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where 5v, SX, Sfi, Se Sp are the functions of the coordinates r] r such that 

\5u\, \5X\, \5fi\ < 1, Se < e , <*P < Po- (13) 

We shall assume that the frame projection of radial velocity of fluid v is the small 
quantity too. 

Let us note that the linearized gauge transformation 

7] = fj + Srj(fj, f), r = r + 5r(fj, r) (14) 



gives 



ds 2 = - a 2 (r/)(l + 5v)drf + a 2 (r/)(l + <5A)dr 2 + a 2 (??)r 2 (l + Sp)dn 2 
= -a 2 {fj) (l + fc/ + 2<ft) + 2%) drf + 2a 2 {fj)(5r + 5rj')dfjdf 2 



+ a 2 (f/) (\ + 5 A + 2<5r' + 2%^) df 2 



+ a 2 (f/)f 2 (1 + d> + 2— + 2-5ri) dQ 2 (15) 
\ raj 

Here the dot and the prime denote d/dfj and <9/<9f respectively, dVL 2 = d6 2 + sin 2 # dip 2 . 
The functions #77 and Sr can be chosen such that 

d(Sr) d(Sri) „ d(<5r) <5r 5u — 5 A .„ ^. 

+ -y^ = o, -y^ - — = — (i6) 

dr) dr dr r 2 

and consequently 5X = Sp,. Omitting the tilde over r], r, 5X and Sfi below we shall 
suppose 

Sfi = 8X. (17) 

Now one can to rewrite the metric in linear (with respect to the perturbations) 
approximation as 

ds 2 = -a 2 { V ) {(1 + Sv)dif + (1 + SX) [dr 2 + r 2 (l + Sfi)dQ 2 } } (18) 

Thus Einstein's equations (jHEJ) are reduced in linear approximation to the following 
four equations 

S\" + 2— - 3-SX + 3°--Sv = -8na 2 Se, (19) 
r a a 2 

Su' + SX' •■ d • d_. f'a d 2 \ 2 Se , nnS 

5X - 2-SX + -Sv + 2 - ) Su = 8ira 2 — , (20) 

r a a \ a a 2 J 7 

8v" + SX" Sv' + 8X' -•• „d • d r . /a d 2 \ . „ 9 Se 

+ 5A-2-5A + -<5t/+ 2 ) <Ji/ = 87ra 2 — , (21) 

2 2r a a \ a ar/ 7 

<JA' - ^z/ = 8vr ^1 + a 2 e ^ (22) 

for the four unknown functions <5A, <fe, v (we take into account 8p = Ss/j). 
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Let us note that as the consequence of the symmetry of problem the perturbations 
described in this paper are scalar ones. The gauge-invariant amplitude of density 
perturbation with the notations of the paper jH] is 

6 3 = 5-3(l + w)i|^)-i^ (23) 

The choice of the isotropic coordinates (|18|) in our paper corresponds to a longitudinal 
gauge. With the notations of the paper [3] this give 

B (o) = H p = ( 2 4) 

By comparing the corresponding expressions in two papers we find 

e g Q {0) = 5Q (0) = Se/e . (25) 
The relations between the other perturbed quantities can be determined analogously 
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From this relations we find that the perturbation quantities introduced in 
this section are actualy physical as well as the gauge-invariant quantities 

e 9 g (0) ,^ 0) g (0) ,$Ag (0) ,^g (0) . 



The difference of equations (|2*U|) (}2"T]) gives 

{SX" + Sv") = - (SX' + Su') . 
r 



(29) 



The solution of this equation is 

5X + Su = d + C 2 r, (30) 
where C\ and Ci are arbitrary functions of the time coordinate 7]. If we require that 

SX = Su = 0, for r — > oo, (31) 
i.e. the space-time is asymptotically homogeneous and isotropic, then 

d = C 2 = (32) 

and 

<5z/ = — <5A. (33) 
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One can to consider the expressions (|19|) and (|22|) as the definitions of 5e and v 
respectively 



3 a / • a 
8irSe = 8ir^5p = — ( 5 X +-5X 



1 



5X" + -SX' 

r 



inv 



a 2 (l + 7)e V a 



(34) 
(35) 



Thus the spherically symmetric perturbations can be described by the metric 

ds 2 = -a 2 { V ) {(1 - 5X)d7] 2 + (1 + SX) [dr 2 + r 2 (l + 5fi)dQ 2 } } , (36) 

where the single unknown function 5X(r],r) is determined by the following combination 
of the expressions (jl9H21j) 

(37) 



5X" + -6X' = j5X + 3 (1 + 7) -5X. 

r a 



4. Particular solutions 

First of all it is necessary to note that the equation (|3Tj) has a solution, which correspond 
to the Newtonian potential [Hj 



5X 



. 2m 
-ov = — , 
ar 



caused by a particle of variable mass 

m( v ) = C 3 a + C 4 a- 3(1+7)/{27) , 

where C3 and C4 are constants. 

If the perfect fluid is a dust, i.e. 

p = (1/7 = 0), 

then the metric perturbations are described by the equation 

d 2 (SX) 7 d(5X) 
da 2 2a da 
The solution of this equation is 



0. 



5X 



,5/2 



+ F 2 



(38) 
(39) 

(40) 
(41) 
(42) 



where F\ and F2 are the arbitrary functions of r. 

The case of the ultra relativistic fluid (e = 3p) has been considered in 7J. The 
corresponding solution has a form 



ar or] [a 



(43) 



where $+(x) $_(x) are the arbitrary functions. 
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5. Boundary conditions and solution near the sound horizon 

The hyperbolic equation (fTTj) describes the outgoing and ingoing waves travelling with 
speed of sound from and into the center of the configuration. Therefore if we like 
to describe a local perturbation of radius r (at the moment 770) we can impose the 
boundary conditions at the sound horizon of this perturbation 



6X\ 



d(SX) 



r— ro— r+ro=0 



Or 



d(SX) 



r—ro—r+To=0 



drj 



0. 



(44) 



r— ro— t+to=0 



where 



r = r?/V7- (45) 

Such boundary conditions are in accordance with the "birth" of a perturbation as a 
result of the redistribution of matter. If we make the substitution 
V(r,r) 



where 



B= 3 - 

H 2 



ra@ ' 

1 + 7 

7 



(46) 
(47) 



and introduce new coordinates 

x = t + r, 

y = T - r - T + r , 

then Eq. (p?Tj) can be rewritten in the form 

B 2 V(x,y) „ / da 



dxdy 



V(x,y)=0. 



(48) 



(49) 



_ adT \T=(x+y+T -r )/2 / 

We shall find the solution of this equation near the boundary y = representing the 
functions (da/(adr)), r x+y+T0 _ rQ y 2 and V(x, y) as a power series 

k 



da 

adT \ T =(x+y+T -r )/2 



E 

k=0 



k\ 







da 



2dr J adr 



\T=(x+T -r )/2 



V(x,y) = Q(y)J2v n (x) y n . 



(50) 



(51) 



n=2 



Here we take into account that the first terms of last series vanish as a consequence of 
the boundary conditions (pHj) . B-function 



e(y) 



0, y < 0, 

1, y>o 



(52) 
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is introduced into expression (|5T|) since the perturbation must be vanish out of the sound 
horizon. Then the equation (JIT?)) is rewritten as 



oo ( 

E 



dV„ „_, „.. A y k+ "+" 



n=2 L 



1 dx 

d \ s da 
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The first two terms of this series give 
dV 2 



0. 



dx 



0. 



dV 3 (3 
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The solutions of these equations are 

^ 2 /3 / da \ 
l/ 3 = ^ 3 -\ — 

6(1-/3) V a ^/| T =(.+ TO -ro)/2 

Here we take into account that 



and 



Thus 



dx 
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(3 / dr 
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da 
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da 



6(1-/3) ^adr 



e(y). 



|r=(x+To-ro)/2 



(x + r -r ) 2 , 

The coefficient 1 / (ror ) in ()46|) must be expand in series of y also 

1 + 



ra^ir 



Thus 



SX(x,y) 
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(53) 



(54) 



(55) 



(56) 



(57) 



(5f 



(59) 



(60) 
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or, if we take into account (|47j) 

(r - r - r + r ) 2 



8X(r,r) 



(r + r - T + r )a^O) 

(9 + 147 + 57 2 ) da(z 



+0 
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- r - t + r ) 
(r + r + r - r Q ) 



r ^2 
^3 + 7-— — ^ 

(r + r - r + r J 
r - r - r + r ) 



2=(r+r+T — r )/2 



For the case ^3 = r 
5X from |H] near the sound horizon 



0(r-r-r o + ro). (61) 
r = 0, 7 = 3, this expression coincides with the expansion of 



6. Conclusion 



We have considered the spherically symmetric perturbations of a perfect fluid on the 
background of spatially flat Friedmann models. It was assumed that the fluid is 
described by the linear equation of state (j3J). It was shown that in longitudinal gauge 
the local perturbations of the energy density, the pressure, the frame projection of the 
radial velocity of fluid (}3"4"1 133j) and the metric coefficients (fTTl l3~3*j) are determined by the 
single function. This function satisfies equation (JHTj) of hyperbolic type and describes 
the outgoing and ingoing waves travelling with speed of sound from and into the center 
of the configuration. The boundary conditions (J4*4"j) on the solution of this equation are 
imposed at the sound horizon of perturbation. The solution is represented as a power 
series (j4*ol loTj) near this horizon. First two terms of series are obtained obviously (J6*T]) . 
The corresponding to these terms solution coincides with the expansion near sound 
horizon early known one in the case of ultrarelativistic equation of state (e = 3p). 
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